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In this paper, we present a high resolution angle resolved photoemission spectroscopy (ARPES) 
study of the electronic properties of graphite. We found that the nature of the low energy excitations 
in graphite is particularly sensitive to inter layer coupling as well as lattice disorder. As a consequence 
of the interlayer coupling, we observed for the first time the splitting of the tt bands by ~0.7 eV near 
the Brillouin zone corner K. At low binding energy, we observed signatures of massless Dirac fermions 
with linear dispersion (as in the case of graphene), coexisting with quasiparticles characterized by 
paraboUc dispersion and finite effective mass. We also report the first ARPES signatures of electron- 
phonon interaction in graphite: a kink in the dispersion and a sudden increase in the scattering rate. 
Moreover, the lattice disorder strongly affects the low energy excitations, giving rise to new localized 
states near the Fermi level. These results provide new insights on the unusual nature of the electronic 
and transport properties of graphite. 
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Carbon, one of the most basic elements in nature, still 
gives a lot of surprises. It is found in many different forms 
( "allotropes" ) - from zero dimensional fullerene, one di- 
mensional carbon nanotubes, two dimensional graphene 
and graphite, to three dimensional diamond - and the 
properties of the various carbon allotropes can vary 
widely^. For example, diamond is the hardest material, 
while graphite is one of the softest; diamond is trans- 
parent to the visible spectrum, while graphite is opaque; 
diamond is an electrical insulator while graphite is a con- 
ductor. What is interesting is that all these different 
properties originate from the same carbon atoms, simply 
with different arrangements of the atomic structure. 

One of the prime examples for the unique proper- 
ties of carbon is graphite. Although graphite is one of 
the most widely studied materials, it continues to pro- 
vide scientists with new and interesting challenges as re- 
cently testified by a wide range of novel discoveries in the 
field, such as novel quantum Hall effect^, room temper- 
ature ferromagnetism.^, metal-insulator-like transition^, 
and superconductivity^. Also, the recent discovery that 
the electronic properties of thin graphite samples can be 
modified by externally applied voltage, in such a way 
that these systems can be switched from n-type to p- 
type carriers^-^, has raised a great potential for carbon- 
based nanoelectronics. The secret to all of these fasci- 
nating phenomena lies in the unique nature of the elec- 
tronic properties of graphite, where the low energy exci- 
tations resemble massless Dirac fermions. Although sev- 
eral transport experiments have shown results in agree- 
ment with the existence of Dirac fermions in graphitei, 
so far direct experimental evidence has been limited. 

In this paper, we report an angle resolved photoemis- 
sion spectroscopy (ARPES) study of the low energy ex- 
citations in graphite and we report signatures for the co- 
existence of massless Dirac fermions and quasiparticles 
with finite effective mass. We investigate the effect of 
interlayer coupling and lattice disorder, which provides 



new insights on the nature of the electronic and tranport 
properties of graphite. 



EXPERIMENTAL TECHNIQUE 

ARPES is one of the most widely used techniques to 
study the electronic structure of solids and is the only 
tool that can directly map the energy-momentum phase 
space of electrons. This is based on the principle of the 
photoelectric effect, explained by A. Einstein in 1905. By 
shining a beam of monochromatic light on the sample, 
electrons ( "photoelectrons" ) are emitted by the photo- 
electric effect (see schematic of a typical ARPES exper- 
iment in Figures l(a,b)) and can be probed as a func- 
tion of their kinetic energy for each given emission angle, 
which uniquely determines the momentum of the pho- 
toelectrons. Therefore, given the conservation laws for 
energy and momentum, the kinetic energy and in-plane 
momentum of the photoelectrons can be related unam- 
biguously to the binding energy and the crystal momen- 
tum of electrons in a solid, providing unique information 
on the direction, the speed and the scattering processes 
of valence electrons. 

For a two dimensional solid, where the electronic struc- 
ture is anisotropic and there is no dispersion along the z 
axis (sample normal direction), the electronic dispersion 
is completely determined by the in-plane momentum ky . 
On the other hand, for a three dimensional sample, tfie 
full momentum information including the out-of-plane 
momentum k^ is important, since the electronic struc- 
ture also depends on k^. Note that in the photoemission 
geometry, ky is a good quantum number while k^ is not. 
Although this makes the interpretation of ARPES data 
for a three dimensional sample more difficult, a simple 
modelling of the photoemission process involving free- 
electron approximation of the final state photoelectrons 
has been successfully applied in ARPES^"^ to obtain full 




FIG. 1: (a) Representation of a photoemission process. Elec- 
trons with binding energy Es are excited to the vacuum and 
their kinetic energy E/c is measured. Here Ef is the Fermi 
energy; E^;ac is the energy of the vacuum level and (j) is the 
work function, (b) Typical geometry of an ARPES experi- 
ment. Photoelectrons are measured as a function of E^ and 
angle G, which can be converted to binding energy Es and 
momentum k for electrons inside the solid. 



momentum information, including the k^ value. In par- 
ticular, at normal emission (k||=0), k^ is related to the 
kinetic energy E/^ by k;2 = ^2m(£'/c -h Vin)/fi^^ where the 
inner potential V^^ is a parameter that can be determined 
from the symmetry of the measured dispersions^'^. 

Under the sudden approximation, ARPES measures 
the single particle spectral function A(k, cj), A(k, oo) = {- 
l/7r)ImG(k, cj), and therefore contains unique informa- 
tion on the self energy I](k, uj) of solids. In the non- 
interacting case, A(k, uo) is simply a delta function (5(k, 
uj). From an experimental point of view, A(k, uj) is di- 
rectly related to the intensity measured in an ARPES 
experiment by I(k, uo) = A(k, a;)f(a;)|M(/c, , where 
f{ij) is the Fermi function, reflecting the fact that pho- 
toemission can probe only the occupied states, and M(k, 
cj), known as matrix element, depends on the electron 
momentum and on the energy and polarization of the in- 
coming photon. Therefore, by measuring the intensity of 
the photoemitted electrons as a function of their kinetic 
energy and emission angles, one can probe the electronic 
structure of the material^. 

A typical two dimensional ARPES spectrum is shown 
in Figure 2. The color scale represents the intensity of 
the emitted photoelectrons, which is plotted as a func- 
tion of their binding energy and in-plane momentum in 
panel a. A traditional and powerful method of analyzing 
ARPES data involves energy distribution curves (EDCs), 
energy scans at constant momentum (panel b). In prin- 
ciple, an EDC analysis can provide direct information on 
the total self energy i;(k, u) of the system in two dimen- 
sions. On the other hand, however, the EDC lineshape is 
complicated due to the presence of elastic and inelastic 
photoelectron scatterings in ARPES and the Fermi func- 
tion, f(c<;). An alternative way to the EDC analysis is to 
analyze the momentum distribution curves (MDCs), mo- 
mentum scans at constant energy (panel c). One of the 
advantages of the MDC analysis is that the lineshape can 
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FIG. 2: (a) Typical ARPES intensity map measured on a 
high temperature superconductor Bi2212, as a function of en- 
ergy and momentum. Different colors represent the inten- 
sity of photoelectrons. The color scale for the intensity maps 
throughout this paper except Figures 5 and 11 is shown in 
the inset of panel a. (b) Energy distribution curves (EDCs) 
for momentum values from ko to ks, indicated by vertical 
tick marks on the top of panel a. (c) Momentum distribution 
curves (MDCs) taken at energies from Ei and Es, indicated 
by horizontal tick marks on the right of panel a. 



be well approximated by a Lorentzian function. However, 
an important caveat of this analysis is that it is often 
based on the unchecked assumptions of linear dispersion 
and momentum-independent self energy. 

In this review, we present high resolution ARPES data 
on graphite taken under different experimental condi- 
tions (polarization and photon energy) with main focus 
on the high symmetry directions. Given its quasi- two di- 
mensional crystallographic structure and the weak Van 
der Waals bonding between graphene planes, graphite 
is an ideal system to be studied by photoemission spec- 
troscopy. High resolution ARPES data have been col- 
lected at beam lines 10.0.1 and 12.0.1 of the Advanced 
Light Source of the Lawrence Berkeley National Labora- 
tory. The energy resolution is from 15 meV to 40 meV, 
and the angular resolution is < 0.3 degree. The samples 
were cleaved in situ in ultra-high vacuum better than 
6.0x10"^^ Torr and measured at a temperature of 25 K. 
The Fermi energy was calibrated by measuring the Fermi 
edge of Au. 

The ARPES data presented in this paper have been 
recorded on highly oriented pyrolytic graphite (HOPG). 
We note that, although in order to perform momentum 
resolved measurements, it is most desirable to have a 
properly aligned single crystal, we have recently shown 
that ARPES study can be performed also on an az- 
imuthally disordered sample^^. Almost paradoxically, 
the complete angular averaging gives a total dominance 
of the dispersions along the two high symmetry directions 
over dispersions in all other directions, and thus surpris- 
ingly sharp dispersions are observed in the ARPES data 
obtained from HOPG samples. More mathematically, 
this has been explained in terms of the van-Hove sin- 
gularities in the azimuthal density of states along the 
high symmetry directions'^, F-M-F' (A-L-A^) and F-K- 
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FIG. 3: (a) Schematic drawing of the crystallographic struc- 
ture of graphite, (b) Three dimensional Brillouin zone (BZ) 
of graphite, (c) LDA dispersion for tt (red) and a (black) 
bands along the high symmetry directions for graphite. The 
TT bands at low energy are the main focus of this paper and 
are highlighted by red color. 



(A-H-L^). For the main interest of this paper, i.e. the 
low energy electronic structure within 1.5 eV below E^?, 
the dispersions are dominated only by those along F-K- 

(A-H-L'), since the dispersions along the other direc- 
tion F-M-F' (A-L-A') lie at much higher binding energy 
(larger than 2.4 eV)ii. 



BAND STRUCTURE AND LOW ENERGY 
EXCITATIONS IN GRAPHITE 

Graphite can be considered as a quasi-2D system from 
both a structural and an electronic point of view. From 
a structural point of view, the basic unit of graphite, 
graphene, is formed by carbon atoms arranged in a hon- 
eycomb lattice. Graphite is formed by infinite layers of 
graphene stacked in an ABAB sequence. In graphite, the 
in-plane C-C bonds are strong due to the small in-plane 
lattice spacing of 1.42 A, while the out-of-plane coupling 
is much weaker due to the large interlayer spacing of 3.37 
A. Thus the interlayer interaction in graphite is domi- 
nated by weak van der Waals force. These structural 
properties impart quasi-two dimensionality to the elec- 
tronic structure of graphite. 

The electronic properties of graphite are highly 
anisotropiGiHJ.2iMiiSii&. In graphite, the strong in-plane 
cr bonds formed by 2s, 2pa^ and 2p^ orbitals result in 3 <j 
bands at high binding energy. The out-of-plane tt bond- 
ing formed by orbitals is much weaker, and the tt 
bands are at much lower binding energy. Among all these 
bands, only the tt bands cross the Fermi energy E^? at the 
corners of the hexagonal Brillouin zone (BZ). These low 
energy tt bands play the most important role in determin- 
ing the electronic and transport properties of graphite 
and will be the focus of this paper. Figure 3(c) shows 
the theoretical band structure of graphite. We note that 
different from a perfect two dimensional graphene, the 
weak interlayer interaction in graphite results in a split- 
ting of the TT bands in the k2=0 plane (K-F-M-K). While 
in the plane of k2=0.5 c* (H-L-A-H), where c* is the re- 



FIG. 4: (a) Linear dispersion near Ef, which is expected for 
graphene. (b) Point-like Fermi surface and cone-like disper- 
sion near BZ corner. 



ciprocal lattice constant along the stacking direction, the 
TT bands are degenerate, and the electronic structure in 
this k^ plane strongly resembles that of graphene. De- 
spite the splitting of the tt bands, the overall dispersion 
of the TT band still bears a strong similarity to that of 
graphene. 

The unique low energy electronic structure of graphene 
and graphite^ is of special interest and believed to be re- 
sponsible for various exotic properties observed^^SsifiiiSSi. 
For most solids, the physics is successfully described in 
terms of the nonrelativistic Schrodinger equation. The 
low energy excitations are quasiparticles with a finite ef- 
fective mass, characterized by a finite Fermi surface and 
finite density of state at E^?. On the other hand, graphene 
and graphite are semi-metals and the tt bands are be- 
lieved to disperse linearly and touch only at one point 
near E^? (Figure 4(a)), forming cone-like dispersions near 
each BZ corner (Figure 4(b)). As a consequence, the 
Fermi surface is characterized by "Fermi points" (red 
dots in Figure 4(b)) located at the six corners of the 
hexagonal BZ, rather than a Fermi surface. In this case, 
the low energy excitations resemble relativistic Dirac 
fermions and are essentially described by the relativistic 
Dirac equation with the effective speed of light replaced 
by the Fermi velocity. The density of states is vanishing 
at the crossing point (known as Dirac point), leading to 
many peculiar phenomena in these materials, e.g. novel 
quantum Hall effec1p2iiSii&, magnetic field driven metal- 
insulator-like transition^, unconventional lifetim e^^'^"*^'^^ , 
and spin-Hall effect in graphene—. 



EXPERIMENTAL RESULTS 

1. BAND STRUCTURE AND INTERLAYER 
COUPLING 

Figure 5(a) shows the measured ARPES intensity map 
as a function of binding energy and in-plane momentum 
over a wide energy range (from E^? to -22 eV). In 
this color scale, red is maximum intensity and blue is 
zero intensity. Despite the azimuthal disorder nature of 
the sample, sharp dispersions from two high symmetry 
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FIG. 5: (a) ARPES intensity map taken at 60 eV photon 
energy (k^ ^ 0.28 c*) on HOPG along the gray hne shown 
in the inset. The color scale for this figure and Figure 11 
is shown in the inset, (b) LDA dispersions along two high 
symmetry directions FKM' (blue lines) and FM (red lines) 
stretched by 20% in energy. 




FIG. 6: (a) Angle- integrated intensity curves taken near nor- 
mal emission measured at different photon energies from 40 to 
140 eV. Filled circles mark the peak positions of the tt bands, 
(b) Extracted peak positions from the angle- integrated inten- 
sity curves as a function of kz . The dotted line is the guide to 
the periodicity of the dispersion. From the symmetry of the 
final states detected, the inner potential is determined, (c) 
LDA band structure of the tt bands at k^^O and k^^O.S c*. 
The energies are stretched by 20%. 

directions FKM' and FM directions can still be clearly 
resolved. The overall dispersion is in good agreement 
with the LDA band structure stretched by 20% in energy 
(shown in panel b for comparison), and in agreement with 
previous measurements24i2Si2&. The contributions from 
different orbitals tt, a can be clearly distinguished. We 
now focus on the low energy tt bands between Ep and 
-11 eV, the only bands that cross Ep. 

In order to investigate the details of the low energy 
electronic structure such as the effect of interlayer cou- 
pling, it is important to obtain the full momentum in- 
formation including the out-of-plane momentum k^. To 
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FIG. 7: (a-b) ARPES intensity map measured near the BZ 
corners at photon energies of 43 (k^ ^ 0.35 c*) and 55 eV (k^ 
^ 0.10 c*) respectively. AB and BB label the antibonding 
and bonding tt bands, (c-d) MDCs at -1.2 eV for data shown 
in panels a and b respectively. 



extract the k^ value, we first determine the inner poten- 
tial Yin from the symmetry of the measured dispersion 
along FA (k||=0). Figure 6(a) shows a few examples of 
angle-integrated EDCs over a wide photon energy range 
from 34 to 155 eV, from which the dispersion is extracted. 
There are two features in panel a, a main peak at higher 
binding energy associated with the tt bands, and a hump 
at lower binding energy associated with the a bands. The 
main peak from tt bands shows oscillation between -7.2 
and -8.4 eV as a function of photon energy, or equiva- 
lently k^. The extracted dispersion (panel b) from the 
peak positions can be described as a periodic oscillation 
riding on a linear slope. The linear slope can be ex- 
plained by angle integration over 8 degrees^^ due to the 
angle average mode used for this data set alone among 
the data shown in this paper. The periodic behavior and 
the symmetry of the dispersion enable us to determine 
the inner potential (17±1 eV) and thus extract the k^ 
values. The periodic dispersion bears a strong similarity 
with previously reported data in the literature^^'^^, but 
covers more than twice the k^ range. We note that, as 
previously reported^^'^^, the dispersion shows a period- 
icity of 2 c* rather than c*. This doubling periodicity 
has been observed^^ and can be understood by combin- 
ing the kz dispersion of tt bands shown in panel c and the 
symmetry of the final states detected using the dipole se- 
lection rules^^. More specifically, due to the nonsymmor- 
phic group in graphite, the final states detected changes 
in a repeated sequence of a — P — j — P — a. We note 
that a and 7 states are resulted from the splitting of the 
TT bands, and thus the measured dispersion at the zone 
center is a reflection of the interlayer coupling. 

We now focus on the effect of the interlayer coupling 
near the BZ corners. Figure 7 shows the tt bands as a 
function of ky near the BZ corner for two different k^ 
values, and 0.5 c*. Panel a shows an ARPES intensity 
map of the tt bands taken near the H point {kz=0.b c*). 
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The data show a A-shaped dispersion crossing Ep near 
the apex (H point), similar to the dispersion of the Dirac 
fermions. The Fermi velocity estimated from the 
slope of the dispersion dE/dk=6 eV-A is yf= 0.91±0.15 
m-s~^, similar to reported values on graphene by trans- 
port measurement3i2*i&. Panel b shows an ARPES in- 
tensity map taken near the K point {'kz=0). In this case, 
one can clearly distinguish two A-shaped dispersions, one 
similar to that in panel a, and the other located at a 
higher binding energy. The intensity variation in the 
bonding and antibonding bands is attributed to the pho- 
toemission dipole matrix element^^. The splitting of the 
TT bands can be confirmed in the MDC shown in panel d, 
where two peaks each from the bonding and antibonding 
TT bands are clearly observed. From the MDC dispersions 
shown and the peak positions in the EDC at the K point 
(not shown), this splitting is estimated to be ^ 0.7 eV. 
This is the first clear demonstration of the splitting of the 
TT bands near E^?, while we note that some data in the 
literature2Si22i^ may now be seen as suggestive of this 
splitting. The splitting of tt bands near the K point but 
not near the H point is in agreement with band structure 
calculationi2ii2iiiiiSii&, which also confirms the validity 
of extracting the kz values. Furthermore, the linear A- 
shaped dispersions shown in panel a, strongly resembling 
those of Dirac fermions, are signatures of Dirac fermions. 



II. LOW ENERGY EXCITATIONS AND DIRAC 
FERMIONS 

We now focus on the nature of the low energy ex- 
citations in graphite. These are of fundamental im- 
portance to understand transport properties as well 
as various exotic behaviors observed in graphene and 
gr aphite^ ^ 18 J-^i.?Q , Recently, tranport studies of graphene 
have shown that the low energy excitations behave like 
Dirac fermions. One interesting finding is that the types 
of charge carriers can be tuned from electron to hole by 
applying an electric (magnetic field), due to the vanish- 
ing density of states near Ep^'^. In graphite, transport 
measurements have suggested two main types of excita- 
tioins (massless Dirac fermions and quasiparticles with 
finite effective mass)^, while direct investigation on the 
properties of these low energy excitations has been miss- 
ing so far. 

Figure 8 shows high resolution ARPES data measured 
along AHL' direction. By following the maximum in- 
tensity in panel a, it is clear that the dispersion shows 
a linear behavior. In panel b, we show the raw MDCs 
at different binding energies. In all the MDCs, one can 
clearly distinguish two peaks, a main peak and a weaker 
one. The dispersion extracted from MDC peak posi- 
tions (open circles in panels a,b) shows a linear behavior, 
strongly resembling the behavior of Dirac fermions. The 
Fermi velocity, or the effective speed of light, is deter- 
mined to be 0.8±0.2xl0~^ m-s~^, also similar to values 
reported from band structure calculation^^ and transport 




FIG. 8: (a) ARPES intensity map taken at 20 eV photon en- 
ergy (k^ ^ 0.46 c*) near H. The blue open circles are the dis- 
persions extracted from MDCs. (b) MDCs at energies from 
-0.2 eV to Ef measured from data shown in panel a. The 
schematic drawing on the right hand side of panel b summa- 
rizes the dispersions near H. 
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FIG. 9: (a) ARPES intensity map near the zone corner K 
measured at 50 eV photon energy (k^ ^ 0.08 c*). Due to 
the different matrix elements caused by different photon en- 
ergies and experimental geometries, the intensity in this case 
is strongly enhanced along FK direction. Open circles are 
the dispersions extracted from MDCs. The inset shows a 
schematic drawing of the dispersion near K. (b) MDCs from 
-0.3 eV to Ei?. The schematic drawing on the right hand side 
of panel b summarizes the dispersions near K. 



measurementsiS. By extrapolating the dispersion, the 
crossing point of these two linear bands (known as the 
Dirac point) lies above Ei?, suggesting that the low en- 
ergy excitations near H point are holes, which is in agree- 
ment with transport measurements^'^^. Moreover, with 
the advantage of ARPES to directly probe the momen- 
tum space, the location of the hole pocket as well as the 
Dirac fermions can be determined directly to be near the 
H point. 

To investigate how the low energy excitations evolve 
as a function of k^, we show in Figure 9 the intensity 
map near the zone corner K (kz=0), along FKM' di- 
rection. In contrast to the linear behavior observed in 
Figure 8 near the H point, the dispersion shows clearly 
a parabolic behavior. The parabolic dispersion can also 
be followed by tracking the peak positions in the MDCs 
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FIG. 10: (a) ARPES intensity map taken at 26 eV photon energy (k^ ^ 0.13 c*). (b) EDO Dispersion and (c) EDO half width 
half maximum extracted from the data, showing that there is an ARPES kink at ~ 180 meV. (d) APRES intensity map taken 
at 20 eV photon energy (kz=0A6 c*). (e,f) A[ and Ej2g phonon modes, which might be involved in the kink observed in panels 
a and c. 



(panel b). The parabolic behavior observed here at kz=0 
clearly suggests that the quasiparticles have a finite effec- 
tive mass. Further analysis to extract the effective mass 
IS m progress'^ . 

To conclude this part, the data shown in Figures 8-9 
suggest that the low energy excitations in graphite are: 
massless Dirac fermions near the H point; quasiparticles 
with finite effective mass near the K point^^. 



III. ELECTRON-PHONON INTERACTION 

Electron-phonon interaction is important in under- 
standing the properties of not only graphite, but also 
carbon nanotubes which can be considered as rolled-up 
graphene sheets. First, electron-phonon interaction is im- 
portant in resolving the long debate on the nature of 
Raman D and G peaks observed in carbons^^. More- 
over, electron-phonon interaction might be responsible 
for the breakdown of the ballistic transport expected for 
quasi-one-dimensional carbon nanotubes^^'^^. Recently, 
electron-phonon interaction has been reported in the 
phonon dispersion of graphite, appearing as two Kohn 
anomalies for the mode at F and the A'^ mode at the 
K point^^'^^. On the other hand, evidence of electron- 
phonon interaction on the electronic structure has been 
missing so far. To study the effect of electron-phonon in- 
teraction in the electronic structure, APRES can provide 
direct information, as has been successfully demonstrated 
in metals and high temperature superconductors^^'^^. 
Here we report the first APRES signature of electron- 
phonon interaction in graphite: a kink in both the dis- 
persion and the scattering rate. 

Figure 10(a) shows an ARPES intensity map near the 
K point. The dispersion, while showing a general linear- 
like behavior, shows a non- trivial structure at ~ 180 
meV. This observation is made clear in panel b where 
the ARPES dispersion extracted from EDC peak posi- 



tions is shown. As shown by the dotted line, the dis- 
persion deviates from a linear behavior at energy ^ - 
180 meV. This type of dispersion change is typical of an 
ARPES kink that arises from the interaction of electron 
with a collective mode, but it also could be the result 
of a genuinely non-linear one electron band dispersion as 
discussed above. In order to distinguish the two possi- 
bilities, it is essential to examine the ARPES line width, 
which reflects the inverse lifetime, or the scattering rate. 
If the energy scale is caused by an ARPES kink, then it 
is expected that the line width shows a sudden increase 
at that energy scale, consistent with an increased inco- 
herence at high energy, while in the one-electron band 
scenario, no such abrupt increase is expected. In panel c, 
we show the EDC width (half width at half maximum on 
the Ef side) as a function of energy, clearly demonstrat- 
ing that the ARPES kink scenario is more appropriate. 
Figure 10(d) shows an ARPES intensity map measured 
near the H point. Near E^?, the dispersion is linear as 
discussed in Figure 8, while a stronger kink in the dis- 
persion and a broadening of the peak above the kink en- 
ergy is clearly observed, suggesting that a stronger kink 
is present near the H point. 



Now that we have shown that there is a clear ARPES 
kink signature in graphite, a natural question follows as 
to its origin. We find that electron-lattice coupling is 
a probable source of the kink observed, since the kink 
energy is similar to the energies of the phonons which 
are responsible for the Kohn anomalies in the phonon 
dispersionSS, « 160 meV for A'^ mode (panel e) and ^ 
200 meV for E2g mode (panel f). Further study to iden- 
tify the boson modes responsible for this kink and inves- 
tigate the momentum dependence of this ARPES kink is 
in progress. 
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FIG. 11: (a) First derivative of raw data taken at 55 eV pho- 
ton energy (k^^ ^ 0.11 c*). The color scale is the same as 
Figure 5. Three nondispersive features are observed at en- 
ergies ei (-2.9 eV), e2 (-4.3 eV) and es (-7.8 eV). (b) EDCs 
at r point (thin line) and k||=0.4 (thick line), showing 
the peaks at ei and es are nondispersive. (c) EDCs for the 
momentum range indicated by a horizontal line from ki to 
ki2, marked in panel a. 



VI. LOCALIZED STATES IN GRAPHITE 

Within the graphene plane or each stacked graphene 
layer in graphite, a varying number of imperfections can 
be found, such as vacancies when lattice sites are unfilled 
indicating a missing atom within a basal plane; stacking 
faults when the ABAB sequence of the graphene planes 
is no longer maintained; and zigzag edge occurring near 
unfilled lattice sites^^'^^. Because of the unique electronic 
structure of Dirac fermions and the negligible density of 
states near the Dirac point, graphite is extremely sensi- 
tive to topological defects which can strongly modify the 
electronic structure and the scattering process of quasi- 
particles, thus resulting in important changes in the elec- 
tronic and transport properties. For example, it has been 
predicted^^ that extended defects as lattice dislocation 
can lead to self-doping effects and presence of localized 
states at the Fermi energy. Self-doping effect results in 
electron or hole pocket at rather than a single point, 
as predicted for an ideal graphene system. Zigzag edge 
on the other hand induces localized states near the Fermi 
energy, resulting in a peak in the local density of states 
near the Fermi energy^^-^^ . Finally the presence of disor- 
der strongly modifies the scattering process, resulting in a 
minimum of the scattering rate at finite binding energy^^ 
instead of Ei? as what is expected for Fermi liquid. 

In Figure 11, we introduce disorder /inhomogeneity fea- 
tures that are pronounced only in the case of HOPG sam- 
ples. Panel a shows the first derivative in energy of an 
ARPES map, from Ep to -11 eV. The first derivative is 
a method which allows to enhance the dispersive features 
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FIG. 12: (a) ARPES intensity map taken at 26 eV photon 
energy (k^^ ^ 0.13 c*). (b) ARPES intensity map measured at 
the same conditions as panel a except on a different spot inside 
the sample. The arrow points to the additional intensity near 
Ei?. (c) EDCs at momentum values labeled as a and f3 in 
panels a and b. An additional peak is observed in EDC taken 
at (3. 



as well as rising or falling edges in the data, and thus is 
particularly useful for detecting nondispersive peaks and 
edges. In panel a we can clearly distinguish three nondis- 
persive features at -2.9 eV, -4.3 eV and -7.8 eV, indicated 
by arrows on the same figure. The nondispersive nature 
of these features can be checked by the EDCs shown in 
panels b and c. These features, appearing as sharp ex- 
tended horizontal lines in panel a suggests the presence 
of nondispersive localized states. The energies of these 
nondispersive features occur at the top and bottom of 
the dispersive band and these features are strongly con- 
nected with the dispersive features associated with the 
TT and a bands. Thus we have associated the presence 
of such nondispersive features to the elastic scattering of 
electrons in either the initial state or the final state by 
inhomogeneity or disorder 

We now discuss the more important effect of disorder in 
the low energy electronic properties, as observed in both 
single crystal graphite samples and HOPG samples. In 
order to investigate this effect, we have performed po- 
sition dependent ARPES study with main focus on the 
near Ei? states. Each spectrum is averaged over ^ 100 
/im, the spot size of the synchrotron beam. We find that 
some of the low energy properties of the electronic struc- 
ture are indeed strongly position dependent and we as- 
sociate them with the presence of disordered states. We 
note that similar behavior has been observed in differ- 
ent kinds of graphite samples measured, suggesting that 
these are important properties associated with graphite. 

In Figure 12, we report typical ARPES intensity maps 
near E^?, taken at 26 eV photon energy in different posi- 
tions on the same sample. While panel a show a parabolic 
TT band similar to the feature discussed in Figure 9, it is 
clear that, as we change position within the same sample, 
we observe an additional weakly-dispersive electron-like 
feature, within 50 meV below Ejr (panel b), coexisting 
with the parabolic tt band. The presence of this addi- 
tional electron-like feature in panel b can be clearly ob- 
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FIG. 13: (a) ARPES intensity map taken at 26 eV photon 
energy (k^ ^ 0.13 c*). Red and blue circles are the dispersions 
extracted from EDCs for the electron-like band and tt bands, 
while the dotted black line is the dispersion extracted from 
MDCs. (b) MDC at Ef. Blue, red and gray dotted lines are 
the peaks used to fit the MDC. (c) MDC at -0.1 eV. (d) EDCs 
taken at momenta from ko to ki5 indicated by vertical tick 
marks on the top of panel a. The EDCs at ko to ks are scaled 
by different factors so that all the peaks can be seen, (e-g) 
EDCs at ki, kr, ki2. 



served in the EDO shown in panel c. Here an additional 
weak peak (pointed to by an arrow in panel c) is clearly 
observed. Note that while the presence of the tt band is 
position independent, observed in all the samples studied 
for all positions measured, this additional electron-like 
feature is strongly position dependent, for each of the 
sample studied. 

Figure 13 shows the detailed analysis of this addi- 
tional band near Ep- The details of this electron-like 
band and the tt band are analyzed from the MDCs and 
EDCs shown in panels b-g. Panel b shows the MDC at 
Ei?. A main peak in the center, associated with the tt 
band, and two side peaks associated with the additional 
band, can be distinguished. As the binding energy in- 
creases, the low energy band is no longer present and the 
peaks in the MDC at -100 meV (panel c) are from the 
main tt band. However, the dispersions of the low energy 
bands are difficult to follow in the MDCs, partly due 
to the highly parabolic bands in this region that makes 
an MDC analysis hard to interpret. On the other hand, 
in the EDCs, this low energy feature shows up clearly 
as a small well-defined peak at ki (panel e), ki2 (panel 
g) and a small hump in ky (panel f) and thus we use 
EDC analysis (panel d) to extract the dispersion. The 
extracted EDC dispersions are overplotted in panel a for 
both the low and high energy bands. By fitting the elec- 
tron pocket, we can directly measure the mass of the 
electrons. This gives a value ~ 0.4 mg, which is much 
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FIG. 14: (a) Raw MDCs at energies from Ef to -0.1 eV. The 
dotted lines are the peaks used to fit the MDC. The solid 
lines are the fit curves, (c-d) Comparison of the central peak 
at different energies, between Ei? and -56 meV (panel c) and 
-56 meV and -100 meV (panel d), after subtraction of the 
side peaks. The MDCs are shifted so that the half maximum 
positions on the left are at the same k point, which allows a 
direct comparison of full width half maximum by comparing 
the half maximum position on the right. Interestingly the 
MDC peak near -56 meV has the smallest width. 



larger than that of electrons and holes as measured in 
transportSi^Si^. Also, by estimating the volume of the 
large electron pocket^^, we obtain an electron concentra- 
tion of ~ 8x10^^ cm~^. This electron concentration is 
again an order of magnitude higher than the value re- 
ported by transport measurement a^i^^ . 

We now discuss the origin of this feature. While it 
seems appealing to associate this electron-like feature 
with the electron pocket predicted by band structure 
at k2=0, we note that its position dependence as well 
as the detailed analysis of the effective mass and carrier 
concentration, clearly suggests a different origin for this 
feature. One possible explanation for this large electron 
pocket is due to defect-induced localized states-^^iMi^, 
e.g. states along a zigzag edge. Additional support comes 
from STM, where a peak in the local density of states at 
an energy (^ -0.03 eV) similar to the weakly dispersing 
electron pocket discussed here^^, is observed near zigzag 
edges. 

To gain more insight into the nature of this low en- 
ergy band, it is important to study MDC full width half 
maximum (FWHM), which gives information about the 
inverse lifetime or imaginary part of the self energy. Fig- 
ure 14(a) shows the quality of the MDC fit from which 
the FWHM is extracted. Within the band width of the 
electron pocket (E^ to -40 meV), the MDCs are fitted 
with 3 peaks. For binding energy larger than the band 
width of the electron pocket (-44 meV to -100 meV), the 
MDCs are fitted with 2 peaks only. Panels b and c show 
the comparison of the main peak from the tt band after 
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FIG. 15: (a) ARPES intensity map taken at 26 eV. (b) MDC 
full width half maximum at low energies between -100 meV 
and Ef- An anomalous upturn around -56 meV can be ob- 
served. The error bar from the fit is included, (c) Imagi- 
nary part of the electron self energy, including both the ef- 
fect of electron-electron interaction and the effect of 0.01% 
disorder^^. 

subtracting the side peaks. Surprisingly, the minimum 
MDC FWHM is at ^ -56 meV (see comparisons in pan- 
els b, c), contrary to the Fermi liquid theory where the 
minimum FWHM is expected to be at Ep. 

Figure 15 shows the extracted MDC FWHM. In panel 
b, from -100 meV to -56 meV, the MDC width decreases, 
while near E^?, a sudden increase of the MDC FWHM is 
observed. This sudden increase is present near the energy 
where the large electron-pocket is observed and therefore 
most likely associated with it. It has been predicted that 
in the presence of disorder-induced localized states, the 
imaginary part of the self energy shows a minimum at 
a finite binding energy rather than at Ei? (panel c]^'^. 
The similarity between panels b and c presents an in- 
triguing possibility that the MDC width analysis is also 
indicative of the fact that the electron pocket stems from 
an impurity band induced by defects. However, as we 
noted earlier, one of the caveats of the MDC analysis is 
the assumption of a linear band dispersion relation. The 



breakdown of this assumption in the current case makes it 
necessary a more sophisticated line shape analysis, which 
is in progress. 



CONCLUSION 

In conclusion, we have presented a detailed high reso- 
lution ARPES study of the electronic band structure of 
graphite. We found that, while the overall band struc- 
ture is in agreement with LDA calculation, the low en- 
ergy excitations present new surprises. First of all, we 
report the first evidence that, as a result of the inter layer 
coupling a splitting of the tt bands occurs far from the 
Fermi level, with an overall splitting of 0.7 eV. Second 
we observed that, as we move closer to the Fermi energy, 
the low energy excitations near the H point are Dirac- 
like massless fermions which coexist with quasiparticles 
with finite effective mass near the K point. In addition, 
we report first ARPES signatures of electron-phonon in- 
teraction in graphite. Finally, we also found that the 
electronic structure is strongly affected by lattice disor- 
der resulting in strong defect-induced localized state near 
the Fermi energy, which strongly contribute to the quasi- 
particle scattering process. 
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